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The structure of a superconducting vortex has been studied theoretically for 
a dirty antiferromagnetic superconductor (AFSC), modelling an AFSC as a 
doped semi-metal with s-wave superconducting pairing and antiferromagnetic 
(dielectric) interaction between electrons (holes). It is also supposed that the 
quasiparticles dispersion law possesses the property of nesting. The distribu- 
tion of the superconducting and magnetic order parameters near the vortex 
core is calculated. It is shown that the antiferromagnetic order, been sup- 
pressed at large distances, is restored around the superconducting flux and the 
vortex core is in fact insulating and antiferromagnetic, in stark contrast to 
the normal metal cores of traditional superconductors. Moreover, our model 
calculations predict that as the temperature decreases the flux region of the 
superconductivity and antiferromagnetism coexistence increases. 

PACS numbers:74.20.De, 74.25. Op, 74.25.Ha 



1. INTRODUCTION 

The influence of the antiferromagnetic (AF) correlation in the super- 
conducting state has been attracting much attention recently. In particular, 
the interrelationship between magnetism and superconductivity (SC) have 
been actively discussed in connection with the problem of the symmetry 
type of the order parameter in high-Tc systems (see, e.g., 1 ). It is natural 
also to ask if the close proximity between antiferromagnetic and supercon- 
ducting phases could have any macroscopic manifestations. At the macro- 
scopic level, the external magnetic field penetrates type-II superconductors 
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via normal state metallic inclusions, or vortices, and can reveal the nature of 
the system's ground state. So that, by weakening the superconductivity of 
high-Tc due to magnetic field, the interplay between antiferromagnetic and 
superconducting ordering can be explored. There are predictions that vor- 
tices in high-Tc cuprates differ from those in conventional superconductors. 
E.g., Zhang 2 has developed a unified theory where antiferromagnetism and 
ci-wave SC in high-Tc superconductors has been constructed based on 50(5) 
symmetry. In this theory the AF and SC order parameters are unified into 
a five dimensional vector (a superspin). The 50(5) theory predicts the ex- 
istence of superconducting vortices with antiferromagnetic cores as a spatial 
class of topological solutions. 3 Recent magnetic neutron diffraction data, 4 
scanning tunneling microscopy, 5 muon spin rotation 6 and nuclear magnetic 
resonance 7 measurements on some high-Tc SC have revealed that the ap- 
plied field which imposes the vortex lattice, also induces antiferromagnetic 
electronic correlations. The new finding is that for optimally doped sam- 
ples (superconducting material shows no static magnetism in either zero or 
weak magnetic field) the vortex state has much stronger tendencies towards 
magnetic order than either the normal or superconducting state. Since mag- 
netism extend beyond the vortex core into the superconducting regions of 
the material, these result means that antiferromagnetic correlation and su- 
perconductivity coexist throughout the bulk of the material, and applied 
field has only effect of recovering some of the magnetic correlation. 3-7 One 
view that this means that the band quasiparticle picture cannot account for 
the AF order of the vortices and the experimental facts provide good support 
of the idea of the existence of 50(5) symmetry in high-Tc cuprates. (see, 
e.g., 3 and references therein). 

Recently it was discovered that competition between superconductivity 
and antiferromagnetism also occur in the new classes of quaternary inter- 
matallic superconductors such as the the borocarbide series RN12B2C and 
RNiBC (here R - rare-earth ion) (see, e.g., 8 and references therein). For 
these new systems, the magnetic energies are dominated by exchange in- 
teractions and comparable to the magnetic superconducting condensation 
energies ensure that the interlay between superocnductivity and antiferro- 
magnetism should be on a comparable energetic footing. 

In this report, we study the possibility for antiferromagnetic ordering 
of the vortex based on the model where the Hamiltonian consists of two 
components, one corresponding to the (Cooper pair) condensation conditions 
of equal and opposite momenta and spins and the other to the magnetic 
interaction condition of parallel spins and finite momentum transfer. A 
mean-field approach is taken to study the magnetic and superconducting 
structures of a flux vortex in such AFSC in the so-called "dirty" limit. We 
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show that the external magnetic field stabilizes the AF order around the 
vortex core in a region where SC have been suppressed, i.e., a coexistence 
phase is induced by the applied field. It was also obtained that the vortex 
core is, in fact, insulating in stark contrast to the standard metal cores of 
low-temperature superconductors. 



Selecting the Hamiltonian of an AFSC , we will use the isotropic model 
of a doped semi-metals with Cooper pairing electrons (holes) within each 
(two) of the bands and dielectric (electron-hole) pairing between bands. 9 In 
the self-consistent field approximation we have 



+ [A 15 (r)^V i + + A 2S (r)<^+ + h.c] + A A (r) J» CT + h.c] h) 



Here ip£, ip a and f^, <£V are the creation and annihilation operators of 
electrons and holes in the first and second bands, respectively, whose dis- 
persion law satisfies the condition of nesting £i j2 (p + Q) = £p- The fraction 
of the Fermi surface that is nested, v , will be treated as a variable. As is 
known, the Fermi surface nesting causes the transition to spin density wave 
state. 9 ' 10 This state can be described by dielectric order parameter, which 
characterizes the gap in quasiparticle's excitations. For the case of an an- 
tiferromagnetic (doubled) cell, we take Q = K/2 , where K is the vector 
of a reciprocal lattice of the crystal. The operator £p is understood to be 
h 2 (p — eA) 2 /2m — Ef; A 1 ( 2 )g(r) are the band superconducting order pa- 
rameters; Ayt(r) is the magnetic order parameter; Ep = Ti 2 p F /2m stands 
for the Fermi energy; and A is the vector potential. As one can see, for the 
model (1) the SC correlations have the spin-singlet type, while the magnetic 
ones have the spin-triplet type. 

The system of Gor'kov equations for electrons in the first band will 
be distinguished from the standard system for functions Gn(l,2) = — < 
T T ipi(l)ipf(2) > and F u (l,2) =< T T ipf (2) > by the presence of the 
additional function 3?2i (1 5 2) = — < T T (p^(l)ipt(2) > which describes the 
AF correlations of electrons from different bands. (Here we used standard 
notifications, i.e., T T is the r-ordered operator, 1 = (ri,ri) and so on.) 
Due to the magnetic interaction between the electrons of different bands, 
the Green's functions of the second band with a reversed sign for the spin 
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The model 
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appear in set of equations. The complete system of equations contains 12 
equations (see Ref.ll for details). In this report the case, when total SC 
order parameter of the system has the s-wave symmetry, i.e., the case 9 A15 = 
A25 = As will be studied. We present here the main equations for such a 
system in so-called dirty limit. 12 



3. Quasiclassical equations for dirty AFSC 

There are three characteristic lengths for the system under consider- 
ation: the electrons mean free path t , the SC correlation length £s ~ 
hvp/Aso and the AF correlation length ^a~^vf/ Aao • Here A$o is the 
superconducting order parameter if the antiferromagnetic pairing is absent, 
while A^o is the antiferromagnetic one when the superconducting pairing 
is absent. The £4 can be considered as the characteristic length of electron 
liquid's magnetic modulation or magnetic stiffness, as well. We will assume 
that the metal is "dirty" for both orders; i.e., the mean free path I and the 
correlation lengths £g an d £,a satisfy the condition £ <C (£s,£a) • Using the 
standard procedure, one can obtain 11 the system of equations, that are the 
generalized Usadel equations for the model (1). 

It is convenient 11 to pass from the band superconducting variables to 
the total ones: A s = A 1S + A 2S , and F = F u + F 22 , G u = G 22 = G 
, Q21 = ^12 = 3\ There are tree type of Green's functions: the functions 
G = G(oj, r) and 3 = r) describe the normal excitations in system of 
electrons with antiferromagnetic correlations, while the anomalous Green's 
function F = F(u, r) describes the Cooper-pair condensate. These functions 
are related by the usual normalization condition 

G 2 (^,r) + |F(u;,r)| 2 + |^,r)| 2 = l (2) 

Magnetic and superconducting properties of the system are described by 
generalized quasiclassical equations that can be written in the form: 11 

- \dU{GUF - FVG) = —r-G - loF (3) 

DV(GV%-%VG) = 2lu%-2^-G (4) 

Where II = V + 2z7rA/$o is the gradient-invariant momentum operator, 
$0 stands for the flux quantum, D = lvp/3 is the diffusion coefficient and 
fiuj = huj n = nT{2n + 1) is Matsubara frequency (T is the temperature, k# 
= 1). 
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These equations should be supplemented by formulas for the order pa- 
rameters and relation between supercurrent and vector potential. The self- 
consistency conditions for the order parameters As and A^ are 

A s HT/T C0 ) + 2ttT £ (As/Hu - F) = (5) 

oj>0 

A A ln(r/r JV o) + 2 7 rr^(A A /^-9) = (6) 

Here Tqo is the superconducting transition temperature when antiferromag- 
netic pairing is absent, while Tjvo is the antiferromagnetic transition tem- 
perature when these is no superconducting pairing. The self-consistency 
condition for the vector potential has the usual form (London electrody- 
namics): 

V x (V x A) =—j s = i—eN(0)DT ^ (FHF + h.c.) (7) 

c c lo>0 

where N(0) is the density of states on the Fermi surface. If an antiferromag- 
netic pairing is vanished the Eqs. (2), (3), (5) and (7) restore those for a 
nonmagnetic superconductor. 

Assuming the two order parameters to be spatially homogeneous, the 
model (1) leads to four possible states, each with a specific solutions (As, A^). 
Namely: (i) the normal state (0,0); (ii) the BCS superconducting state 
(As, 0); (iii) a purely antiferromagnetic state (0, Aa); and (iv) a coexistence 
phase where both As and A^ are non-zero. The phase with the lowest free 
energy is stable. The phase diagram for the system can be calculated in 
terms of the relative interaction strengths of the two orders (the relationship 
Tco/Tno) and the fraction of the Fermi surface that is nested (the results 
will be published elsewhere. 13 ) We will not analyze here all the phases but 
only the BCS one when in the sample bulk the antiferromagnetic order is 
suppressed. That means that the boundary conditions for Eqs. (2)- (4) 

Go = hu/y/(hu) 2 + A| , F = A so /-\/(M 2 + A|o, = (8) 

match the solutions to the Green's function at large distances. Here A50 is 
the superconducting order parameter if antiferromagnetic pairing is absent. 

4. Basic equations in the vortex geometry 

To study the vortex it is convenient to work in the gauge in which F 
and As , and 9 and A^ are real. The coordinate dependence of all functions 
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is then only through the distance r from the axis of the vortex and we may 
define ip = ip(u n ,r) = ip n and 9 = 9(u n ,r) = 9 n by such a way that the 
Eq.(2) is fulfilled: 

G(lo, r) = cos ip cos 9, F(u, r) = sin p cos 9, Ss(oj, r) = sin 9 (9) 

Due to the cylinder symmetry of a vortex the lines of A are circles around it 
and the relations VA = 0, AV(/? = and AV# = took place. By making 
use these relations, the equations of motions for ip and 9 can be obtained 
from Eqs. (3) and (4) in the form: 

£ 2 {cosfl(vV) - 2sin0(Vy>)(V0)} = 

1 9 /2ttA\ 2 „ A s (r) hoj , . 

^ sin2<pcos0 + — cos p - —— sm p (10) 



£ 2 {cos p(V 2 9) + sin 9 [sin <p cos cos (p cos 9{S/p) — 2 sin <p sin 0(V<p)(V#)]} 

hu . A A (r) , . 

sinS cos p cos 9 (11) 



A50 A 50 

where £ = (D/2Aso) 1//2 stands for the SC correlation length. The equation 
for the vector potentials now is took the form 

V x (V x A) = --QA, (12) 

where the diamagnetic response function Q reads 

Q = A7re 2 N{0)DT/hc x 2vr ^ sin 2 p cos 2 9 (13) 

o;>0 

At large distances from the flux axis ( r — > 00 ), in accordance with the 
boundary conditions (8) we have ip — > (^oo = arct(/(A5o/w) and 6* — > = 0. 
Then Q — > Qo = c/47rA 2 where A is the London penetration depth is given 
by the usual dirty limit formula 14 

A- 2 = ge 2 iV(0)Dtanh^ (14) 

In the limit of extreme type-II SC the solution to Eqs. (10), (11) and (12) 
may be developed by a series of powers of Ginzburg-Landau parameter 1/k 2 
(k = A/£) . We will calculate the results up to 0(1/k 4 ) terms. Note, that 
for k 3> 1 and £ < r the expression for vector potential is dominated by the 
largest term, which, as follows from Eq. (12), is 

he 

Mr) = -^K^r/x). ( 15 ) 
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Here K\(x) is the imaginary Bessel function of first order. Note that if 
Aa = 0, 6 n = our formulas restore those considered by Watts- Tobin and 
Waterworth. 15 

For the purposes of computation, it is convenient to choose R = r/£ as 
the independent variable and to obtain 

d 2 <p n 1 dip n 1 _ 2 zsii-oi \ A s (i?)cos(^ n ftu; n sin</? n 
— o + — — TTT = -K SinZWooAi (R/k) : — h-7 tt - (16) 

dR 2 R dR 2 ^ 11 M A 50 cos n A so cos n v 7 

d 2 6 n ld6 n ( A S (R) . ^hco n A A (R) 

~im + = l ~a sm ^« + -t— cos ¥>„} sin 6> n cosfl„ (17) 

dR z R dR A S o A S o A S o 

Here in the terms of oc kT 2 , the zero-order values of A(r), <p(r) and 9{r) are 
used. Accordingly, the self-consistency conditions now are determined by 

A S (R) = 2vrT sin Vn cos 6 n {ln ^— + 2vrT £ l/^}" 1 (18) 

n>0 C0 n>0 

A A (i?) = 2vrT ^ sin #„{ln + 2vrT ]T l/fe^}" 1 , (19) 

n>0 JjV0 n>0 

and A50 is determined by the relation (18) with cos9 n = 1. 

We solved these set of equations numerically by an iteration procedure. 
The results of our calculations are presented in Figs. 1-3, where the SC and 
AF order parameters near the vortex core, Ag(R) and Aa(R) , respectively, 
are shown for different values of the Tco/Tno ratio, for a given fraction of 
the Fermi surface that is nested and at different temperatures. The solid line 
represents a coordinate dependence of the superconducting order parame- 
ter As(R); the dashed line depicts the antiferromagnetic order parameter 
Aa(R) behavior. As it is seen in the figures, the size of the "antiferromag- 
netic" vortex core is on the order of the superconducting coherence length. 
I.e., the model predicts that the vortex core is in fact insulating and anti- 
ferromagnetic. Moreover, as T/Tco decreases the region near the flux core 
where superconductivity and antiferromagnetism to coexist increases. 



5. Summary 

We investigated the structure of an isolated single vortex for a dirty 
metal with competing superconducting (a spin singlet pairing) and antifer- 
romagnetic (a spin triplet pairing) interactions. We shown that the external 
magnetic field stabilized the antiferromagnetic order around the supercon- 
ducting core, in the region where superconductivity has been suppressed. It 
is obtained that the vortex core is insulating and antiferromagnetic, in stark 
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contrast to the standard normal metal cores of traditional superconductors. 
The main our nontrivial result is that as T/Tco decreases the region where 
superconductivity and antiferromagnetism to coexist increases. 

In recent years, one can see an explosion in the field of vortex dynamics, 
especially in the materials of high-Tc SC. As is well known, inertial mass 
of a fluxion play an important role in vortex dynamics. There are various 
mechanisms contributing to the flux inertial mass, and the important one 
is due to variation of local electronic density of state in Abrikosov vortices. 
The metallic core with states very close to the Fermi energy predicted by 
Caroly et. a/. 16 is the standard picture for discussion of inertial mass of 
the vortex. However, the static, long-ranged, magnetism associated with 
the vortex core may have nontrivial influence upon transport properties of 
magnetic superconductors in an applied magnetic field and should be taken 
into account for systems with competition between superconductivity and 
antiferromagnetism. 

The authors are grateful to V. Chabanenko and N. Hayashi for stimu- 
lating discussions. 
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Fig. 1 Radial coordinate dependence of the vortex superconducting 
order parameter Ag(R) (solid lines) and antiferromagnetic order parameter 
Aa(R) (dashed lines) for k = 10 and different values of T/Tco = 0.05, 0.5, 
and 0.8. T co = T m ; v < 1. Fig. 2 Same as in Fig. 1 for T co = 0-5T m , 
v < 0.5. Fig. 3 Same as in Fig. 1 for T co = 2T m , v < 1. 



